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HI 21cm absorption at z ~ 2.347 towards PKS B0438-436 
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ABSTRACT 

We report the detection of redshifted HI 21cm absorption in the z ~ 2.347 damped Lyman- 
a absorber (DLA) towards PKS B0438-436, with the Green Bank Telescope. This is the 
second-highest redshift at which 21cm absorption has been detected in a DLA. The absorp- 
tion extends over ~ 60 km s _1 and has two distinct components, at z = 2.347477(12) 
and z — 2.347869(20). A similar velocity structure is seen in optical metal lines, al- 
though the peak absorption here is offset by ~ 11 km s _1 from the peak in the 21cm 
line. We obtain a high spin temperature T s ~ (886 ± 248) x (//0.58) K, using a cover- 
ing factor estimated from 2.3 GHz VLBI data. However, the current data cannot rule out 
a low spin temperature. The non-detection of CO 3-2 absorption places the upper limit 
N co < 3.8 x 10 15 x (T x /10) cm" 2 on the CO column density. 
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1 INTRODUCTION 

The highest HI column density absorption systems in quasar spec- 
tra, the damped Lyman-a absorbers (DLAs), are important in the 
context of galaxy evolution as the y are beli eved to be the precursors 
of present-day galaxies IWolfe et alll9 86). Understanding the na- 
ture of a typical DLA and its evolution with redshift is perhaps the 
key to understanding normal galaxy evolution. However, despite a 
tremendous observational effort over the last two decades, physical 
conditions in high z DLAs are still the subjects of much debate (see 
IWolfeet all 1 2005) for a recent review). 

HI 21cm absorption studies of DLAs towards radio-loud back- 
ground sources provide crucial information on the nature of the ab- 
sorbing galaxy (e.g. Kanekar & Briggs 2004). 21cm lines are usu- 
ally optically thin and trace the kinematics of the neutral gas; they 
can thus be used to obtain the velocity field of the absorber and its 
HI and dynamical masses iBrig gs et alJl200ll) . In some cases, the 
gas kinetic temperature can be directly measured with high res- 
olution 21cm spectroscopy (e.g. lLane et alJl2000h . Comparisons 
between the redshifts of 21cm and metal-line or molecular ab- 
sorption allow estimates of changes in fundamental constants (e.g. 
iTzanavaris et all2005l) . Finally, the 21cm equivalent width can be 
combined with the HI column density to obtain the spin tempera- 
ture T s of the absorbing gas (e.g. Kan ekar & Chengaluj2003h : this 
contains information on the distribution of the HI between phases 
at different temperatu res, one of the issues of con troversy referred 
to above iKanekar & Chengaluj2003LlWolfe et all2003l) . 

A problem with using 21cm absorption to probe physical 
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conditions in high z DLAs is that there are, at present, only 
three absorbers in t he literature with confirmed detections of ab- 
sorption at z > 1 IWolfe & Dayi3ll979l rWolfe & Briggslll98lt 
Wolfe et al. 1985), with one further tentative detection at z ~ 3.38 
(Briggs et a U I19971) . Despite deep searche s, the vast majority of 
system s only yield upper limits on T s (e.g. IKanekar & ChengaluJ 
2003). This is very different from the situation at low red- 
shifts, z < 0.8, where 14 intervening 21cm absorbers are 
known (e.g?lBrown & Robertsll 19731 IKanekar & ChengalujEoOll 
Lane & Briggs 2002). We have therefore been carrying out deep 
searches for 21cm absorption in a large sample of high redshift 
DLAs with the Green Bank Telescope (GBT) and report here the 
detection of a new 21cm absorption system in the z ~ 2.347 DLA 
towards the z ~ 2.863 QSO, PKS 0438-436 lEllison et alfcOOll) . 



2 OBSERVATIONS AND DATA ANALYSIS 

The initial GBT observations of PKS B0438— 436 were carried 
out in lanuary 2004, with the PF1-450 MHz receiver. These used 
the GBT Spectral Processor as the backend, with two linear po- 
larizations and a 5 MHz bandwidth centred at 424.5 MHz and 
sub-divided into 1024 channels (giving a velocity resolution of 
~ 6.9 km s _1 , after Hanning smoothing). A standard "position- 
switched" mode was used, with ten-minute On-Off cycles made 
up of ten-second integrations. System temperatures were measured 
during the observations using a noise diode. The total on-source 
time was 30 minutes. 

These observations resulted in the tentative detection of an 
absorption feature at ~ 424.32 MHz, close to the expected red- 
shifted line frequency. We hence re-observed PKS B0438— 436 in 
August 2005 to confirm the absorption, at higher spectral resolu- 
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Figure 1. Final GBT HI 21cm spectrum towards PKS B0438— 436 from August 2005, with flux density (in mJy) plotted against barycentric velocity (in km/s), 
relative to z = 2.34736. [A]: The spectrum at the original resolution of 1.73 km s _1 . [B]: The spectrum (dashed line) and 2-gaussian fit (solid line) after 
smoothing by 5 channels (and resampling) to a resolution of ~ 8.6 km s -1 . The QSO continuum flux density has been subtracted from both spectra. 



tion. The observing and data recording procedures were the same 
as those of January 2004, except that a bandwidth of 1.25 MHz, 
centred at 424.3 MHz was used, sub-divided into 1024 channels 
(i.e. a spectral resolution of ~ 1.73 km s _1 , after Hanning smooth- 
ing). The total on-source time was 80 minutes. 

The data were analysed in the AIPS++ single dish package 
DISH, using standard procedures. No radio frequency interfer- 
ence (RFI) was found at or near the expected line frequency. Af- 
ter initial data-editing, the spectra were calibrated and averaged 
together, to measure the flux density of the background quasar, 
using a telescope gain of 2 K/Jy. This gave flux densities of 
7.5 ± 0.7 Jy in January 2004 and 7.1 ± 0.7 Jy in August 2005, 
at 424 MHz, where the errors include those from con fusing sources 
in the primary beam; (for comparison. lLarge et alj (1981) obtain 
8.12 ± 0.25 Jy at 408 MHz). A second-order spectral baseline was 
then fit to each 10-second spectrum (during the process of calibra- 
tion) and subtracted out; the frequency range covered in the fit was 
~ 423.7 - 424.7 MHz, for both data sets. The residual 10-second 
spectra were then averaged together, to produce a spectrum for each 
epoch. The August 2005 spectrum was found to not be entirely flat 
at the end of this procedure; a second order polynomial was hence 
fit to this spectrum, excluding line channels, and subtracted out to 
produce the final spectrum. This has a higher sensitivity, a finer 
resolution and a better spectral baseline than the spectrum of Jan- 
uary 2004, and will therefore be used as the final 21cm spectrum 
towards PKS B0438— 436 in the discussion below. 

We also carried out a search for redshifted CO 3-2 absorption 
from the z ~ 2.347 DLA with the 3-mm receivers of the Australia 
Telescope Compact Array (ATCA), in October 2005. The anten- 
nas were in the EW214 configuration, with a longest baseline of 
214 m. A bandwidth of 64 MHz was used for the observations, 
centred at 103.315 GHz and divided into 128 channels; this gave 
dual linear polarization spectra with an effective velocity resolution 
of 1.77 km s~ . The observations consisted of 10-minute scans on 
PKS B0438— 436, interleaved with 2-minute scans on the bright 



nearby source PKS B0454— 463. The telescope pointing and sys- 
tem temperature calibrations were updated every hour; we obtained 
above-atmosphere T sys in the range 300 — 350 K in the different an- 
tennas, at the observing frequency. The planet Mars was observed 
at the start of the run to calibrate the absolute flux density scale. 
The total on-source time was 210 minutes. 

The ATCA data were analysed in the software package 
MIRIAD, using standard procedures. The complex gains and 
bandpass shapes of the ATCA antennas were derived from the vis- 
ibility data on PKS B0454— 463 while the flux density scale was 
bootstrapped using the visibility amplitudes on Mars measured on 
the shortest (30-m) baseline (Mars is resolved on the longer ATCA 
baselines). The channel velocities towards PKS B0438— 436 were 
then shifted to the barycentric frame and the calibrated visibilities 
averaged to derive the final spectra. The source flux density was 
measured to be 0.36 Jy, at a frequency of 103 GHz. 



3 SPECTRA AND RESULTS 

The final redshifted HI 21cm spectrum towards 
PKS B0438-436 from August 2005 is shown in Figs. [flA] 
and [B], with flux density (in mJy) plotted against barycentric 
velocity (in km/s), relative to z = 2.34736, the red shift of the 
strongest metal line absorption lAkerman et aljE005h . We adopt 
this as the systemic redshift of the DLA, indicated by the dashed 
vertical lines in the two panels of the figure. Fig.QA] shows the 
HI spectrum at the original velocity resolution of ~ 1.73 km s _1 , 
with an RMS noise of ~ 10.4 mJy per channel, while Fig.QFS] 
shows the spectrum smoothed by 5 channels (and resampled), at a 
final resolution of ~ 8.6 km s _1 , with an RMS noise of ~ 4.7 mJy 
per resampled channel. 

The stronger 21cm component of Fig.0B] is clearly detected 
in both observing runs, with the correct Doppler shift due to the 
Earth's motion; there is no doubt about its reality. While the weaker 
21cm component (at ~ +45 km s _1 ) was only detected in the 
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Figure 2. ATCA 13 km s" 1 resolution CO 3-2 spectrum towards 
PKS B0438— 436, with flux density (in mJy) plotted against barycentric ve- 
locity (in km s — 1 ), relative to z = 2.34736. The QSO flux density has been 
subtracted from the spectrum. The dashed vertical line shows the CO 3-2 
frequency corresponding to the peak 21cm line redshift. 

higher sensitivity spectrum of August 2005, the data showed no ev- 
idence for RFI and the feature is detected at > 5a level. Moreover, 
as we will discuss later, both the redshift and the velocity width of 
this feature are in excellent agreement with those of a component 
detected in the Fell A2374 transition. We conclude that the second 
21cm component is also likely to be real. 

The 21cm absorption towards PKS B0438— 436 extends, be- 
tween nulls, over ~ 60 km s" 1 , with an equivalent width of 
~ (0.216 ± 0.027) km s" 1 . The solid line in Fig.[HB] shows 
a two-component fit to the spectrum, with the components mod- 
elled as gaussians. The peak redshifts of the components are z\ = 
2.347477(12) and z 2 = 2.347869(20) and their peak optical 
depths (the ratio of line depth to continuum flux density, using the 
low optical depth limit; PKS B0438-436 is unresolved by the GBT 
beam) are n = (6.2 ± 0.8) x 10~ 3 and r 2 = (3.3 ± 0.7) x 10" 3 
(including errors in the source flux density). They have fairly large 
velocity spreads (FWHMs of ~ (23.3±2.7) km s _1 and ~ (18.2± 
4.5) km s -1 , respectively), suggesting either a blend of many cold 
narrow components or that bulk/turbulent motions are important. 
Alternatively, if the absorption arises in a warm phase, with the ve- 
locity spread primarily due to thermal motions, the implied kinetic 
temperatures of the two components are < (11865 ± 2750) K and 
< (7240 ± 3580) K. 

Fig.|2|shows the final ATCA redshifted CO 3-2 spectrum to- 
wards PKS B0438— 436, with flux density, in mJy, plotted against 
barycentric velocity, in km s _1 , relative to z = 2.34736; the 
dashed line indicates the peak of the 21cm absorption. The spec- 
trum has been smoothed to (and resampled at) a velocity resolution 
of ~ 7.25 km s -1 and shows no evidence for absorption, with an 
RMS noise of ~ 8.6 mJy per 7.25 km s" 1 . The 3u upper limit on 
the optical depth is r < 0.058, per ~ 10.2 km s -1 . Curiously, 
however, a weak emission feature can be seen in the spectrum at 
the redshifted CO 3-2 frequency. While this has only ~ 3<j sig- 
nificance in a single channel (after smoothing by 3 channels and 



resampling), it is somewhat tantalizing, especially as it is present 
(albeit at even weaker levels) in the two independent polarizations 
and occurs at the 21cm redshift. We discuss the implications of a 
positive detection in the next section. 



4 DISCUSSION 

For optically thin 21cm absorption, the HI column density Nm (in 
cm -2 ), 21cm optical depth T21 and spin temperature T s (in K) are 
related by the equation 

TVhi = 1.823 x 10 18 [T S //] j r 21 dV , (1) 

where the profile integral is over velocity (in km s _1 ) and the cov- 
ering factor / gives the fraction of the radio source that is ob- 
scured by the foreground cloud. The detection of 21cm absorp- 
tion in a DLA, where the HI column density is known from the 
Lyman-a line, thus allows an estimate of the spin temperature 
of the absorbing gas (assuming that the HI column density mea- 
sured towards the optical QSO is the same as t hat towards the ra- 
dio source). As has often bee n emphasized (e.g. lCarilli et all 19961 
IChenealur & KanekarlEoQcl) . the spin temperature thus obtained 
is the column-density-weighted harmonic mean of the spin tem- 
peratures of different HI phases along the line of sight. In the 
case of the z ~ 2.347 DLA, the HI colum n density is iVjai = 
(6 ± 1.5) x 10 20 cm" 2 lEllison et aj200ll) . The measured 21cm 
equivalent width of / r 2 idV = (0.216±0.027) km s" 1 then yields 
T s = [(1527 ± 428) x /] K, where the error is dominated by the 
uncertainty in TVhi . 

In order to estimate the covering factor, it is important to first 
determine whether both 21cm absorption components of Fig.QFS] 
arise against the quasar core. Figs.[3]A] and [B] show a comparison 
between the HI 21cm profile (histogram) and the Fell A2374 and 
Znll A2062 profiles of lAkerman et alj fcOOSh . The Fell and 21cm 
profiles are similar, with Fell absorption seen at the velocities of 
both 21cm components; this suggests that the latter indeed arise 
from the line of sight towards the optical nucleus. 

The covering factor / can thus be assumed to be the frac- 
tion of source flux density contained in the core and is best es- 
timated from VLBI observations of the radio source at the red- 
shifted 21cm line frequency; unfortunately, no such observations 
of PKS B0438— 436 exist in the literature. The lowest frequency 
with published VLBI observations is 2.3 GHz, where most of the 
source flux densi ty stems from two compact components, separated 
by 35 mas iPreston et alll989h . Of these, the northwest component 
is more compact than the southeast one (~ 5 mas v/s ~ 18 mas) 
and is likely to be the nuclear core, especially since it is also very 
compact in higher resoluti on (< 2 mas) 4.8 GHz VLBI and VSOP 
images IShen et aHll998tlTingav et alJl2002h . The core has a flux 
density of ~ 2.5 Jy in th e 2.3 GHz VLBI im age, ~ 58% of the 
total 2.3 GHz flux density iPreston etalJl989l) . Unfortunately, the 
source is highly variable and it is thus not possible to determine 
its spectral index or its core flux density at 424.3 MHz from the 
2.3 and 4.8 GHz VLBI data, as these observations were carried out 
at very different epochs, many years before our observations. For 
now, we are forced to estimate the covering factor solely from the 
2.3 GHz VLBI information. Assuming that only the core is com- 
pletely covered gives / = 0.58 from the 2.3 GHz data, and, hence, 
a high spin temperature, T s = (886 ± 248) K. Of course, it is 
quite possible that both 2.3 GHz source components are covered 
by the DLA, since their separation of 35 mas corresponds to a lin- 
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Figure 3. A comparison between the GBT HI 21cm profile of Fig.QB] (histogram) and the [A] Fell A2374 and [B] Znll A2062 profiles (solid points), with 
all spectra in arbitrary units. The redshifts of the two 21cm components are shown by the two dashed vertical lines while the arrows at the top of each panel 
(marked A, B, C and D) indicate the 4 Fell components of Table 1. See text for discussion. 
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Table 1. The Fell and Znll column densities and [Zn/Fe] for the four 
strongest Fell components (in log arithmic units), with the [Zn/Fe] values 
relative to the solar abundances of Lodders 12003). Components C and D 
have velocities in good agreement with those of the 21cm components. 



ear scale of only ~ 290 pc at z = 2.347 (using an LCDM cos- 
mology, with fi m = 0.3, Qa = 0.7 and Ho = 70 km/s Mpc^ 1 ). 
In such a situation, with / ~ 1, the spin temperature would be 
even higher, T s ~ (1527 ± 428) K. Conversely, it is possible that 
the spin temperature is significantly lower than this, if the core has 
an inverted spectrum and is the only source component covered 
by the DLA. We note that the other confirmed z > 1 HI 21cm 
absorbers, in the DLAs towards QSO 1331+170, PKS 1157+014 
and PKS 0458— 020, also have relatively hig h spin temperatures, 
T s > 500 K (e.g. lKanekar & Chengaluj20ol . 

The average metallicity of the z ~ 2. 347 DLA is qui te high 
for these redshifts, [Zn/H] = -0.68 ± 0.15 jAkerman et all2005l) . 
Of course, the lack of velocity information in the Lyman-a line 
precludes measurements of the HI column densities (and hence, 
the metallicities) of individual absorption components. We have 
used VPFIT to fit the Znll and Fell transitions detected in the VLT- 
UVES spectrum of A kerman et alj 120051) : this has a resolution of 
~ 7 km s~ x and S/N ~ 25 near the Znll and Fell wavelengths. 
Table 1 lists the velocities (relative to z — 2.34736), redshifts, b- 
parameters, Fell and Znll column densities and the dust depletion 
(as measured by [Zn/Fe]) of the four strongest components in the 
DLA, obtained from the fits. The first three components (A - C) in 
the table contribute to the strongest optical absorption (peaking at 



v ~ km s _1 in Fig. [3] ~ 11 km s _1 blueward of the stronger 
21cm absorption), while the fourth is associated with the weaker 
21cm absorption, at v ~ 45 km s _1 . The weaker 21cm compo- 
nent is not detected in the Znll lines; 3<j upper limits on its Znll 
column density and dust depletion are listed in Table 1. Although 
the optical absorption has been decomposed into sub-components 
during the fitting procedure, the modest S/N of the UVES spectrum 
and the complex kinematics of the metal lines introduce a consid- 
erable amount of degeneracy. Errors on the column densities and 
redshifts of individual components are therefore sizable. However, 
the total column densities are well constrained, with our values in 
excellent agreement with those of Akerman et al. (2005), despite 
the very different fits. We hence combine the three optical com- 
ponents at v ~ km s" 1 to obtain [Zn/Fe] = (0.62 ± 0.14), in- 
dicating a fairly high dust depletion compared to other DLAs at 
this redshift. Conversely, the component at v ~ 45 km/s has a 
much lower implied dust-to-gas ratio, [Zn/Fe] < 0.26. Fig. 8[A] of 
IWolfe et alJ <2005l) shows that all DLAs with [Zn/Fe]< 0.4 have 
[Zn/H]< —1. Although this result is obtained from [Zn/H] and 
[Zn/Fe] values averaged over the entire profile (i.e. not from in- 
dividual components), it suggests that the v ~ 45 km s _1 gas has 
a low metallicity. The v ~ km s _1 gas should then have a fairly 
high metallicity, to account for the high average metallicity of the 
absorber. The DLA thus appears to contain some gas that is rich 
in dust and metals and with a relatively large 21cm optical depth, 
as well as a second component that is poorer in heavier elements 
and with weaker 21cm absorption. The very different depletions 
seen in the two components imply that the absorber has a rather 
non-uniform interstellar medium, as recently seen in other high z 
DLAs iDessauges-Zavadskv et ai]|2006l) . In fact, the Fell A2374 
and Znll A2062 lines at v ~ km s _1 have quite different shapes, 
suggesting that even the sub-components within this complex have 
significantly different physical conditions. 

Comparisons between 21cm and optical redshifts in a statis- 
tically large absorber sample can be used to probe the evolution 
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of fundamental constants IWolfe et alJll976h . The most sensitive 
result is that o flTzanavaris et alT]200^ . who compare the redshifts 
of strongest 21cm and metal-line absorption in 8 DLAs to constrain 
changes in the quantity x = g p /j,a 2 , where a is the fine structure 
constant, g p , the proton gyro-magnetic ratio and [i = m e /m p , the 
electron-proton mass ratio. However, it is by no means essential 
that the strongest 21cm and metal absorption both ar ise from the 
same absorption component; for example, IChengalur & KanekaJ 
(2000) noted that the deepest UV absorption in the z ~ 2.04 DLA 
towards PKS 0458—020 is offset from the strongest 21cm compo- 
nent but in good alignment with the secondary 21cm component. 
Similarly, both panels of Fig.|3|show that the strongest optical ab- 
sorption in PKS B0438-436 (at z ~ 2.347360) is offset from the 
strong est 21cm a bsorption by ~ 11 km s" 1 . In the approach of 
iTzanavaris et alJ l2005h . these offsets would imply either evolution 
in the fundamental constants or an intrinsic velocity offset between 
the 21cm and metal lines. However, it can be seen from Fig. |3 A] 
and Table 1 that there are weaker Fell components (C and D) in 
PKS B0438— 436 much closer to the redshifts of both 21cm com- 
ponents. While component C (at v = +9.5 km s _1 ) is severely 
blended with the complex at v ~ km s _1 , making it difficult to 
determine an accurate redshift from the present spectrum, the fit to 
component D at v ~ +45 km s _1 is quite stable, yielding the red- 
shift z\ = 2.347864(5) (note that the RMS-error of the wavelength 
solution is ~ 250 m s _1 , a factor of two smaller than the error 
from the fit). Comparing this to the redshift of the nearest (weaker) 
21cm component, Z2 = 2.347869(20) (from the 2-gaussian fit), 
we obtain [Ax/x] = Az/(1 + z) = (-0.15 ± 0.62) x 10" 5 , 
where Az = 21—22 and z is the mean of the redshifts; this is 
consistent with the null result of no evolution. Note that this er- 
ror does not include possible systematic velocity offsets between 
the 21cm and UV redshifts, which can only be addressed with 
statistically large samples. However, it is clear that the results ob- 
tained from such comparisons between different species (and using 
lines at very different wavelengths) critically depend on the pre- 
cise details of the comparison. The assumption that the strongest 
absorption in the two species arises in the same component in- 
troduces an extra source of error, that could well dominate the 
~ 6 km s - 1 "intrinsic" offset be tween 21cm and UV redshifts 
obtained by iTzanavaris et all 120051) . leading to a bias in favour of 
a detection (especially in small samples). It may thus be better to 
compare the redshifts of the nearest simple (i.e. unblended) metal 
and 21cm components, after using high resolution spectroscopy to 
decompose the profiles into their components, as has been done 
here. While this might have the opposite bias, towards a non- 
detection of evolution, it could be considered a test of the null hy- 
pothesis. In any case, care must be taken in the interpretation of the 
results from such comparisons. 

The ATCA non-detection of CO 3-2 absorption in Fig. [2] 
places limits on the CO column density in the z ~ 2.347 DLA. 
Assuming a velocity width of 10 km s _1 gives the 3cr limit Nco < 
3.8 x 10 15 x (T x / 10) cm , similar to earlier limits in DLAs (e.g. 
ICurran et alj20 04f). A conversion factor of 10 5 from Nco to Nh 2 
iLiszt & LucasfcOOOh gives Nh 2 < 3.8 x 10 20 x (T x /10) cm" 2 , 
i.e. a molecular fraction / < 0.6. 

The possible CO 3-2 emission feature seen in Fig. [2] if real, 
corresponds to a molecular mass of Mh 2 ~ 4.2 x 10 10 M©, an 
exceedingly large value (using the Galactic conversion factor of 
4.6 M(T)(K km/s pc 2 )" 1 from the CO line luminosity L' co to 
mass; Solo mon et aflll997l) . CO emission has never before been 
detected from a DLA. We note that the velocity width of the feature 
is quite narrow, ~ 30 km s _1 ; if real, this would be surprising for 



such a large molecular mass and would require the galaxy to be 
quite close to face-on to the line of sight. 

In summary, we have used the GBT to detect HI 21cm absorp- 
tion in the 2 ~ 2.347 DLA towards PKS B0438-436. We obtain a 
high spin temperature T s = (886 + 248) x (//0.58) K, estimating 
the covering factor / from 2.3 GHz VLBI observations. 
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